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ABSTRACT 


A two  server  loss  system  is  considered  with  N classes  of  Poisson 
arrivals,  where  the  service  distribution  function  and  server  prefer- 
ences are  arrival  class  dependent.  The  stationary  state  probabilitie 
are  derived  and  found  to  be  independent  of  the  form  of  the  service 
distributions. 
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by 

Ronald  W.  Wolff  and  Charles  W.  Wrlghtson 


Introduction 

Emergency  service  systems  have  been  modeled  and  analyzed  as  systems 
with  multiple  servers  in  parallel  for  which  no  queue  is  allowed  [1,3]. 

The  classical  model  for  this  type  of  system  is  the  Erlang  loss  model  which 
assumes  Poisson  arrivals  and  general  independent  service,  where  arrivals 
finding  all  servers  busy  depart  without  receiving  service.  For  this 
model  it  is  well-known  that  the  limiting  (stationary)  state  probabilities 
satisfy  Erlang's  formula  independent  of  the  form  of  the  service  distribu- 
tion function  [6]. 

Carter,  Chaiken  and  Ignall  [1]  considered  the  problem  of  designing 
response  areas  for  two  emergency  service  units  with  the  dual  objectives  of 
minimizing  the  average  response  time  to  calls  for  service  and  equalizing 
the  workload  between  the  two  units. 

Chaiken  and  Ignall  [2]  derived  the  theoretical  results  for  the  par- 
ticular queueing  model  used  in  [1]  to  determine  the  response  areas  f >*•  the 
emergency  units.  Specifically,  Chaiken  and  Ignall  assumed  an  M /C/2  model 
with  no  queue  allowed  (a  two-server  Erlang  loss  system)  vhere  arrivals 
have  preference  for  particular  servers.  There  are  two  arrival  classes, 
j c {1,2}  , where  the  classes  form  independent  Poisson  streams  at  rates 
A^  and  , respectively.  Class  j arrivals  choose  server  j if  the 
system  is  empty.  Service  times  are  independent  and  identically  distri- 
buted (the  same  distribution  for  both  classes  of  arrivals).  Under  these 
conditions  they  derived  expressions  for  the  stationary  probabi 1 ities 


and  PQ1  , where 
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Pjq  ■ Prob  (server  1 is  busy,  server  2 Is  free), 

Pqj  * Prob  (server  1 Is  free,  server  2 is  busy),  and 

Pjq+P  «P  ■ Probability  that  one  server  is  free. 

This  paper  generalizes  their  lesults  as  follows:  (1)  an  arbitrary 

number  of  arrival  classes  is  considered;  (2)  each  arrival  class  is  allowed 
to  have  its  own  service  distribution  function;  and  (3)  server  preferences, 
which  still  depend  on  the  arrival  class,  are  expressed  probabilistically. 

In  addition,  the  method  of  proof  is  simpler. 

In  reference  to  the  original  problem  of  designing  response  areas  for 
emergency  service  units,  it  is  reasonable  to  expect  that  calls  for  emergency 
service  from  different  sections  of  an  urban  area  will  have  different  dis- 
tributions for  the  length  of  time  required  to  service  the  calls.  Thus,  this 
generalization  of  tee  above  model  would  be  expected  to  extend  the  applica- 
bility of  the  basic  approach  and  to  improve  the  accuracy  of  results. 

Basic  Mod e 1 and  Notation 

The  basic  model  is  that  of  the  M/G/2  queue  with  no  queue  allowed.  That 
is,  there  are  two  servers,  identically  distributed  exponential  interarrival 
times,  and  identically  distributed  service  times,  where  these  random 
variables  are  mutually  independent. 

Enumerate  arrival  classes  1,2,  ...,  N where  for  class  j arrivals 
we  define: 


X = 


<1,  = 


the  arrival  rate  of  all  j^bs,  i.e.,  the  reciprocal  of  the 

mean  interarrival  time,  where  0 < X < <*>  , 

the  probability  an  arrival  belongs  to  class  j (assumed 

independent  of  the  interarrival  times  and  the  classification 

N 
\ 
t 

i- 


of 


other  jobs),  ]>  q = 1 

j=l  J 
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aJ 


(i.J) 


(t) 


j , the  rate  of  j-arrivals, 
the  service  distribution  of  a j-arrival, 
the  service  time  of  a j-arrival,  P(Sj  < t)  * Cj(t)  , 

‘j  E(V  • 


N 

£ q.G  , the  service  distribution, 
j-1  J J 

the  service  time,  i.e.,  a random  variable  with  distribution 


G , 

N 

l 0 ,«  X E(S)  , 
j-1  J 

the  probability  that  a j-arrival  chooses  server  ] if  the 
system  is  empty  upon  arrival, 

the  probability  that  a j-arrival  chooses  server  2 if  the 
system  is  empty  upon  arrival, 

state  of  system,  i-0  means  server  1 is  free,  i-l  means 
server  1 is  busy,  j-0  means  server  2 is  free,  j-1 
means  server  2 is  busy, 

Prob  (system  is  in  state  (i,j)  at  time  t)  , 
stationary  probability  of  system  being  in  state  (i,j)  = 
lim  P„  (t)  . 


If  exactly  one  server  is  available  upon  arrival,  the  arrival  is  served 
by  that  server  regardless  of  the  arrival  class.  Arrivals  finding  both 
servers  busy  depart  immediately  and  are  not  served. 


Main  Results 

The  main  results  in  this  paper  are  described  in  the  following 


Theorem ■ 

For  the  model  described  in  the  previous  section. 
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!>00  = 1/(1  + c + . “/ 2)  , 


(1) 


poi  * poo  H (1“j)*j  +r2/2)  /(i  + ) ' 


p = p “12 
11  00  ' ’ 


independent  of  the  form  of  the  service  distributions. 

To  prove  the  theorem,  observe  that  if  arrival  preferences  and  the 
distinction  between  servers  are  ignored,  the  system  can  be  viewed  as  an 
Erlang  loss  system  with  two  servers,  Poisson  arrivals  at  rate  1 , and 
identically  distributed  service  times  with  distribution  function  G . 
For  this  system  it  is  well-known  that  the  stationary  probabilities  for 
the  number  of  busy  servers  satisfy  Erlang's  formula  [6]: 


I>  = 


lim  P0 ( t ) 
t— 


• 1/(1  + r + c2/2)  , 


P,  = 


P„  = 


lim  Pj  (t ) 

t -*cc 


lim  ( t ) 
t-*® 


= , / ( 1 + c + r / 2 ) , 


2 2(1  + c + . 2 12)  , 


where  P^ ( t ) = Prob  (j  servers  busy  at  time  t)  , j =0,1,2  . Thus, 
P00  = P0  ’ P11  = P2  ’ anc^  t*,e  Pr°blen'  *s  reduced  to  finding  P^  and 

P01  ' where  P10  + P01  ■ P1  • 

For  the  above  system  with  different  classes  of  arrivals  and  server 


preferences,  one  quantity  of  interest  is  the  proportion  of  time  each  server 
is  busy.  This  quantity  is  dependent  upon  the  arrival  rates,  service  dis- 
tributions and  server  preferences  of  the  arrival  classes.  Clearly, 
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P 


si 


+ P 
+ P 


11 

11 


where  Pgj  - the  proportion  of  time  server  J As  busy,  J * 1,2. 

The  approach  used  in  this  paper  to  derive  P^q  and  is  based 

oa  two  well-known  facts  in  queueing  theory.  First,  the  proportion  of 
time  that  a server  is  busy  serving  a given  stream  of  arrivals  is  equal 
to  the  product  of  the  arrival  rate  and  the  expected  service  time  for  that 
arrival  stream.  This  can  be  viewed  as  a corollary  of  L * XW  [4], 

Second,  the  stationary  probabilities  have  a time  average  interpretation: 

* the  proportion  of  time  the  system  is  in  state  (i,j)  in 
steady-state, 


and,  in  steady-state,  the  system  behavior  found  by  Poisson  arrivals  is 
identical  with  the  time  average  behavior,  i.e.,  the  proportion  of  Poisson 
arrivals  that  find  the  system  in  state  (i,j)  is  equal  to  the  proportion 
of  time  the  system  spends  in  state  (i,j)  [5,6,7].  Similar  approaches 
have  been  employed  to  derive  other  quantities  of  interest  in  various 
types  of  queues  with  Poisson  arrivals  [5,7]. 

The  proportion  of  time  server  1 is  busy  can  be  divided  into  the  propor- 
tion of  time  server  1 is  busy  serving  each  of  the  different  classes  of 
arrivals  [4].  For  class  j , the  proportion  of  time  server  1 is  busy 
equals  the  product  of:  (1)  the  rate  of  class  j arrivals  served  by  server 

1,  X . say,  and  (2)  the  expected  class  j service  time.  Thu*, 

J * 1 


si 


N 

l 

j-1 


K(S4) 
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where  >. . . is  simply  times  the  proportion  of  j-arrivals  that  are 

J • 1 J 

served  by  server  1. 

Therefore , 


,i  ■ £ ypoi  * »)poo>  E<y 


cP01  * r00  J,  °jcJ 


We  also  know  that 


si 


- P10  + P 


11 


01 


+ p 


11 


Solving  for  , we  obtain 


P01 


p + p - P 
1 11  00 

1 + c 


l Yj 
1-1  J J 


Substitution  of  the  values  for  Pq^qq)  * * an£*  P,>^Pll^  * rom  (1) 

yields  PQ1  in  the  theorem.  From  P^q  a P^  - P^  , we  obtain  P ^ . 
This  concludes  the  proof. 


Consideration  of  Various  Extensions 

For  a two  server  loss  system  with  Poisson  arrivals,  arrival  classes 
with  different  server  preferences,  and  service  time  distributions  which 
depend  only  on  the  arrival  class  (arrival-dependent  service  times), 
we  obtained  simple  results  which  are  independent  of  the  form  of  the  arrival 
class  service  distributions.  Can  Erlang's  formula  be  extended  further? 
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Certain  obvious  extensions  were  investigated: 

(a)  three  or  more  servers, 

(b)  server-dependent  service  times,  and 

(c)  preference-dependent  service  times,  e.g.,  in  the  two  server 
pure  preference  (non-probabil ist ic)  case,  the  service  time 
distribution  depends  on  whether  the  arrival  is  served  at  the 
server  of  first  choice. 

For  (a),  Erlang's  formula  will  still  hold  for  the  distribution  of 
the  number  of  busy  servers.  Thus,  the  same  approach  is  feasible; 
however,  there  are  insufficient  equations  to  solve  for  all  of  the  state 
probabilities.  In  addition,  it  wa  verified  that  the  state  probabilities 
depend  on  the  form  of  the  arrival  class  service  distributions.  For 
three  servers,  numerical  solutions  for  the  state  probabilities  under 
exponential  service  were  found  to  be  different  from  those  under  hyper- 
exponential service. 

Erlang's  formula  for  the  distribution  of  the  number  of  busy  servers 
does  not  hold  for  cases  (b)  and  (c).  Also,  it  was  verified  that  :.he  state 
probabilities  are  distribution  dependent  (again,  by  comparison  of  ex- 
ponential and  hyperexponential  service  distributions). 

In  conclusion,  none  of  the  extensions  considered  in  this  section  are 
valid,  i.e.,  results  do  depend  on  the  form  of  service  distributions.  This 
is  unfortunate  because  of  numerous  applications  for  these  models  and  other 
variations  on  loss  systems.  Exact  results  will  almost  certainly  be  hard 
to  obtain  and  will  dej end  on  service  distributions  in  complicated  ways. 
However,  in  some  cases,  partial  results  do  hold  and  are  easily  obtained, 


■Mfc 


t li 


8 


e.g.,  for  case  (a).  Furthermore,  (limited)  numerical  experience  indicates 
that  dependence  on  the  form  of  service  distributions  is  small.  If  this 
is  true,  numerical  results  under  exponential  service  may  be  good  approxi- 
mations. Of  course,  it  would  be  desirable  to  be  able  to  quantify  this 
last  assertion. 
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